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Abstract We consider a molecular single electron transistor coupled to a vibrational 
mode. For some values of the bias and gate voltage transport is possible only by ab- 
sorption of one ore more phonons. The system acts then as a cooler for the mechanical 
mode at the condition that the electron temperature is lower than the phonon tem- 
perature. The final effective temperature of the vibrational mode depends strongly on 
the bias conditions and can be lower or higher of the reservoir in contact with the 
oscillator. We discuss the efficiency of this method, in particular we find that there is 
an optimal value for the electron-phonon coupling that maximizes cooling. 
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1 Introduction 

Cooling the vibrational mode of nano-mechanical oscillator is a major challenge on 
the road to the observation of quantum effects in mechanical systems. Recently several 
proposals for cooling mechanical modes have been formulated [Tll2 ,3,4,5,6,7,8, 9l llOI 
I11II12II13U14U15U16) and many have been reahzed experimentally [I7a8.19.20.21) . A 
particularly efficient method is the so called sideband cooling (see for instance Refs. O 
|6]). The general idea of side-band cooling is to exploit a detuning between the frequency 
of a forcing field (for instance a laser source) and the resonant frequency of a cavity 
mode. The cavity mode is coupled to a mechanical oscillator of resonant frequency ui. 
When the detuning is negative (the source-field frequency is lower than the cavity mode 
frequency) and close to uj the increase of the occupation of the cavity mode can be 
realized by the absorption of a phonon of the oscillator. This process removes energy 
from the mechanical mode and thus cools the oscillator. For positive detuning one 
observes the opposite phenomenon of heating due to the coherent pumping of energy 
in the mechanical mode. 
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Fig. 1 An example of the system considered in the text, a suspended nanotube contacted to 
two leads and coupled capacitively to a gate. 

In this paper we consider a similar phenomenon that can take place in a typical 
nano-electromechanical device: a molecular single electron transistor coupled to a me- 
chanical oscillator (see Fig. [l]). This system has been widely studied in the literature 
[221I23II24II25II26II27II28II29II30II31) . In particular it has been shown that the coupling to 
the mechanical mode can induce a suppression of the current at low bias both in the 
molecular case '25 , 29 , 30] (when the size of the central island is so small that a single 
electronic level participates in the transport) and in the metallic classical case (when 
the level spacing in the island is negligible with respect to the bias voltage) |32) . The 
suppression of the current is known as Frank-Condon blockade or phonon-blockade, and 
it has a very simple classical origin: When the island is charged the force acting on it 
can displace the island itself leading to a different gate voltage (Vg) seen by the device. 
This corresponds to a shift of the Coulomb diamonds in the Vg-V plane {V being the 
bias voltage) by a small amount in the Vg direction. The current can then flow only in 
the regions corresponding to the intersection of the sequential transport regions of the 
shifted and un-shifted diamonds. By a simple geometrical construction one finds that 
at low voltage current is suppressed |32) . In the quantum limit huj ~ fc^T, eV, (where 
T is the environment temperature, kg is the Boltzmann constant and e is the electron 
charge) the shift of the diamonds is discrete and each step corresponds to the absorp- 
tion or the emission of phonons. The resulting Coulomb diamonds are characterized 
by emission [251129] or absorption sidebands [33] ■ Very recently the emission sidebands 
predicted by these theories have been observed in suspended carbon nanotubes |34j . 

In this paper we discuss under which conditions it is possible to exploit the phonon- 
assisted tunneling to cool the oscillator. A necessary condition is that the electronic 
temperature is lower than the temperature of the environment coupled to the oscil- 
lator. This fact can be realized either because the oscillator is coupled better than 
the electrons to a hot environment, or by a preliminary cooling of the electrons by 
existing techniques of electrons cooling |35ir7]. In this paper we explain how one could 
sequentially cool also the mechanical mode by exploiting the properties of the single 
electron transistor to optimize cooling. We find that it exist regions in the Vg-V plane 
where the cooling efi^ect is maximum, these are the regions just outside the conducting 
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diamonds within a distance in energy of the order of hui. We also find that in order to 
achieve a good cooUng efficiency the electron-phonon coupfing constant must not be 
too large, since then the same mechanism at the base of the phonon blockade blocks 
transport at low bias thus reducing also the number of absorbed phonons. 

The paper is organised as follows: In Section [2] we describe in details the system 
and the theoretical approach in terms of a master equation for the phonon distribution. 
In Section[3]we give an analytical description of cooling for small bias voltage and near 
the border of the conducting region in the Vg-V plane. In Section |4] we discuss the 
numerical solution of the master equation for several values of the parameters. Section 
Ogives our conclusions. 

2 System model 

We consider large molecules (like carbon nanotubes) or quantum dots connected to two 
leads and in presence of a gate voltage (see Fig. [l]). If the Coulomb charging energy 
is sufficiently high this system acts as a single electron transistor, where transport 
can be finely tuned by varying the gate voltage. We assume that the presence of an 
excess electron in the central island generates a force on a nano-oscillator, that can 
be a vibrational mode of the molecule itself or a nearby mobile and gated lead. We 
consider the regime where the oscillation frequency is larger or of the same order of 
the temperature, the bias and the gate voltage. This regime can be easily realized for 
small molecules, where the oscillating frequency can attain the THz [221126] . in larger 
molecules and for metallic oscillators this can be more difficult since for instance a long 
carbon nanotube the oscillation frequency can be in the MHz range |36II37) rendering 
the system classical [32] . A minimal model for such a system is the Anderson-Holdstein 
Hamiltonian that has been considered by many authors [38II39II25II27II40] : 

k,a k,a 

Here a is the lead index (L or R), c, c^, d, and S are the electron annihilation and 
creation operators for the leads (c) and local orbital (d). e^.^^ = ^fc — /ia is the dispersion 
relation of the single-electron bands in the lead a with jj,a being the chemical potential. 
The index k labels the single electron states in the two leads, eg is the energy eigenvalue 
of the local orbital and can be shifted by the external gate voltage (cfr. Fig. [T}. The 
operators a and annihilate and create phonons, A is the electron-phonon dimension- 
less coupling. We consider the model for spinless electrons for simplicity. Inclusion of 
spin along with large onsite Coulomb repulsion leads to qualitatively similar results, 
as far as the Kondo temperature ^ fiWjfc^Te, with Te the electronic temperature. 
We also assume a wide band limit, so that the tunneling induced level width is simply 
Fa = Tvuata, where Ua is the density of states and ta is the tunneling amplitude in 
lead a. The exact definition of the electron phonon coupling A depends on the explicit 
geometry and material parameters of the device. Values larger than unit are possible, 
they correspond to a displacement of the mechanical oscillator over a distance larger 
than the zero point motion length when an electron is added in the dot. Very recently 
A has been evaluated theoretically for the elongation modes of a suspended nanotube 
[34] . In good agreement with the observed magnitude of the Frank-Condon sideband 
the authors found that A « 1.3, thus exactly in the intermediate range relevant for the 
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present paper. The mechanical state for fc^Te, -Ta ^ Tiui has been discussed in Ref. (30l 
141] for the molecular case and in Ref. [42ll43ll32ll44] for the metallic case. (The case of 
a single tunnel junction is discussed in Ref. [45].) The frequent hopping of electrons 
in and out the central island are a strong source of decoherence for the oscillator and 
the reduced density matrix is diagonal in the position and momentum basis of the 
oscillator itself. A classical approach is thus appropriate in this case. Here we consider 
the opposite limit of fast phonons: fiuj ^ min(_r^,_rfl) [25ll27|[29ll46) . The mechanical 
oscillator on average performs many revolutions between two electron hopping events. 
The oscillator reduced density matrix in first approximation is then diagonal in the 
base of the energy eigenvectors (phonons). The time evolution of the system can then 
be described by a (Pauli) master equation for the probability Pnp that the central 
island is occupied by n electrons (with n = or 1 in our case) and the oscillator is in 
the eigenstate of the free oscillator Hamiltonian with eigenvalue nTiu). The rates from 
one state to the other can be calculated with Fermi golden rule in second order in the 
tunneling coupling and all orders in the interaction constant A. The standard approach 
is to perform a Lang-Firsov canonical transformation to eliminate the coupling term 
proportional to A from the Hamiltonian in favor of a renormalised tunneling amplitude 
ta — > tae~^^°'~'^ and renormalised energy level eo — ^ = ^0 ~ X^huj. The calculation 
of the transition matrix elements is then straightforward. One obtains (see for instance 
[2511291147]): 

= ^\Mpp' I'fie'o + {p - p)nLo - fic) (2) 

Wp^p,^a = ^Wpp- ? [l - /(^O - {P - P)f^ - Ma)] , (3) 

where Wp~^p/ ^ is the transition rate from the occupation state n to n' with simultane- 
ous change of the phonon state from ptop'. In Eqs. JH) and @ /(e) = 1/(1 + e'^/'''^^= ) 
is the Fermi distribution of the electrons in the metallic leads and Mppi is the matrix 
element for a phonon transition. The latter is given by the overlap integral of two 
harmonic oscillator wavefunctions spatially displaced by Ax — Xy^2h/muj, with m the 
oscillator mass: 




Mp,p,| = A-^e--/^^L^''(n• (4) 



The absolute value of the matrix element depends only on p = min(pi,p2) and P = 
max(pi,p2)- Here Lp are the generalized Laguerre polynomials. We neglect cotunneling 
events, the reader can find a detailed discussion of their contribution to the current in 
Ref. IS]- 

The mechanical oscillator is coupled to an external environment at a given temper- 
ature Tm that we will assume in general different from the electronic temperature Te 
entering the Fermi distributions. If the environment correlation time is shorter than all 
the other time scales of the problems {P/h and following Ref. |29| we can model 
the effect of the environment with the time scale r necessary to reach the phonon 
equilibrium state. The master equation takes then the form: 



— -PnpWnp + ^ Pn'p'Wpi^p a - - j Pnp - Pp'' Pnp' ] (5) 
n'p' a \ p' 
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where Wnp = E„'p'a W;^^^[^ and Pp"« = ^-vn^lkBT^(x „ ^-^u^jknT^^ 
ical equilibrium distribution. The stationary state of the system can now be obtained 
by solving Eq. ((5} for dPnp/dt — 0. Electronic transport can drive the oscillator out of 
equilibrium and in general the resulting effect can be either to increase or to decrease 
the average energy of the oscillator. In the next Sections we discuss the properties of 
the stationary solutions of Eq. © and the consequences for the oscillator state. 

3 Cooling by phonon absorption 

In this section we provide an analytical solution for the stationary state of the device 
and discuss how the bias condition on the single electron transistor can influence the 
thermal state of the mechanical device inducing either cooling or heating. 

Since the stationary solution obtained from Eq. (O is not in general an equilibrium 
distribution strictly speaking we cannot define a temperature for this state. Never- 
theless one can always define the average phonon energy Eav = ^'^X/npP^"P- 
particular, for an equilibrium distribution one has Eav = fiuj/{e^'^/^^'^'^ — 1), thus an 
effective temperature can always be defined by inverting it: 

keT^ff =hu/\n{l + huj/Eav) ■ (6) 

In the following to describe the thermal state of the oscillator we will discuss either the 
average phonon energy or the effective temperature of the oscillator defined by Eq. 

Let us now discuss the form of the solution qualitatively, and in particular the 
contribution of the different tunneling processes. For simplicity let us assume that 
the electronic temperature kgTe <^ hui ~ k^Tm, then the Fermi functions in ([2)| 
and ([3]) are, within a good approximation, step functions. We define the bias voltage 
eV ~ ~ I^R s-nd the gate voltage eVg — {^j^ + A'fl)/2 — Eq in terms of the chemical 
potentials of the two leads. This notation allows to recover the standard description for 
the single-electron transistor as a function of V and Vg. For simplicity in the following 
we concentrate on the V > Q and Vg > Q case. (The results are symmetric in both 
variables if the tunnel junctions are symmetric.) In this case the electrons are on average 
transmitted from the left reservoir to the right reservoir. Sequential tunneling from the 
left to the right reservoir is possible only if some of the M^^ZTp' l ^^'^ ^p^p' R ^''^ 
non-vanishing at the same time. For Te = the relevant conditions are: 

fVKp^Vi/O if eVg^eV/2<-r^Ap 

\ %^::p',_R / if -eVg + eV/2 > huoAp ' ^ ' 

In Fig.[2]we show the different regions of non- vanishing rates, their intersections divide 
the plane into identical rhombi. For Ap = (p' — p) = Eq. ((?]) defines the standard 
Coulomb diamond conducting region in the Vg-V plane, delimited by the two lines 
V — ±2Vg (dashed in Fig. [5]). The conditions Q for Ap ^ divide the plane in regions 
where the value of the rate is constant. The variations are limited to the borders of 
these regions and are discontinuous. 

To draw a simple picture of the mechanism of cooling we consider in the following 
the four relevant cases where the V and Vg are in the regions indicated in the Fig. [5] 
by A, B, C and D . 

Region A is the region of maximal suppression of the current due to the electron- 
phonon interaction |25II29) . It is characterized also by a bistability [321141] since in 
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Fig. 2 (Color online) Regions of non-vanishing rates in the Vg-V plane. The shaded regions 
with W^^" indicate the regions where the rate W^^^,^ 7^ with Ap = p — p' . The dashed 
line indicates the border of the sequential tunneling transport regime for the non-interacting 
case A = 0. 



the limit of large A both the occupied and unoccupied states are blocked and thus 
stationary; the final state of the device depends then on the previous history. In order 
to study the stationary state we begin by assuming that only the probabilities P„o 
and Pni are significantly different from zero. This assumption has to be verified at 
the end of the calculation self-consistently and is equivalent to demanding that the 
stationary state of the oscillator is essentially the ground state with a small occupation 
of the first excited state. We can then solve the problem with only four states: _P = 
{PqO) -Poll ^lOi ^ll}- According to the conditions Q only few rates are non vanishing 
in region A, specifically; 



w' 



-+1 

■*Q.,R ' 



Wtp 



-^0 



-+0 



Wc 



(8) 
(9) 



where Wc = roX^e~^ /h is the relevant cooling rate (/q = Fj^ = Pr)- Note a crucial 
point, all the rates increasing the number of phonons (and thus leading to heating) are 
vanishing. Only the elastic transport rates (|8}, or the cooling rates (|9]l are present, thus 
one can easily anticipate that the system will act as a cooler in this regime. The master 
equation can be written in matrix form as follows dP_/dt = MP_ where the matrix M 
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reads 



f-WTO-Pr/r PqVt Wto 2Wc 

PI'^/t -Wti-2Wc - PqVt Wt 



Wtq 2Wc -Wto-P^Vt P^yr 

\ Wti PiVr -WTi-2Wc-PnVT, 

(10) 

and we have introduced the terms due to the coupling of the oscillator to the external 
environment. The stationary solution is given by the eigenvector of vanishing eigenvalue 
of the matrix M that reads 

P = {Po"« + 2tWc, Pq' + '^'^'^c, Pi"^} ■ (11) 
It gives for the occupation of the excited state Pq = Poq + Piq: 

peq 

This expression has two simple limits: For r ^ the environment of the oscillator 
dominates and the system remains in thermal equilibrium with his own bath. For 
T — > oo, the system reaches its ground state, since Pi — > 0. In general the relevant 
parameter to determine which of the two cases is realized is the product tW^. It 
is maximum for A = 1 that thus constitutes in this region the optimal value of the 
electron-phonon coupling for cooling efficiency. Specifically, for this value of A, Pi 
reads P^'^/(l + 27or/e?i). The condition for efficient cooling is then simply /qt/?!, ^ 1. 

To understand better the behavior of the device let us now consider the three other 
regions next to the phonon blockade area: regions B, C and D, as shown in Fig. |2l In 
region B the situation is very similar to region A, the non vanishing rates are 

W°^Il = W^^Ir = Wtp (13) 

W?^Il = W?^lu = Wi^ro.L = W/i^o%, = Wc. ( 14) 

Again we have no heating terms, the only difference is in the elastic transport terms 
(|13l) . We find that the stationary solution for Pi has the form (|lll) as in the the region 
A. The final effective temperature of region B thus coincides with that of region A. 

Let us now consider region C, this is an important case since one can show that 
actuaUy aU the blocks delimited by the two lines eV = 2eVg and eV — 2eVg + 2h,uj for 
eVg > —2eV + TiLu (indicated as C in Fig. ^ have exactly the same behavior. The list 
of non- vanishing rates reads: 

W^ZIl = <rp,fl. = Wtp (15) 
W^ZIl = W^ZIr = 1^11^0%, = Wc (16) 
W^ZIl = Wh = Wc. (17) 

The situation is different from case A and B, since now one phonon-emission rate is 
non vanishing [Eq. (|17p ]. Still there are three phonon absorption terms (|16p and we can 
expect that the final result will be not far from the ground state. Also in this case it is 
possible to find the stationary solution analytically, but its form is rather cumbersome, 
we prefer to present here only the result for large r: 

^^/^" - "^W^ - (V^ro)e^^2 + A^)/2A^ (18) 
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We thus find again that for large r the single electron transistor behaves as a good 
cooler for the oscillator. It is important to notice that this result holds for large eV, 
since all the blocks along the threshold of current behave in the same way. This holds 
for any value of A, but the optimal range for cooling is A « 1. We can now check the 
self-consistency of the solution. The runcation of the space to the first excited state is 
justified by the fact that Pi ~ /o/r vanishes for tFq oo. 

Let us finally check what happens when the system is biased in region D, where 
the sequential transport channel is open. Only the elastic rates differ from the set of 
non- vanishing rates of region C [Eq.: (|15|l is substituted by the following Equation] 

<rp,L = <i:p,fl, = . (19) 

This difference is very important, since now it is possible to let an electron out of the 
central island by elastic tunneling (W^g^g /?) and let an other electron in by emitting 
a phonon (Wq^^ i). This process heats up the device and is in competition with the 
cooling processes. In region C the only way for an electron to exit the central island is 
by absorption of one phonon by the oscillator (Wj^ZTo _r), thus also if an electron could 
enter the device by emitting one phonon VKqZTi L the net balance of the two processes 
is always of no energy exchanged. Since all the other processes imply cooling of the 
oscillator the stationary solution is a cold state. In region D this is no more true, and 
we can verify it by finding the stationary state for t = oo. We find that Pi does not 
vanish in this limit: 

p ^ Wr^jWc + Wti) , . 

1 4WtoWti + W"c(3W"to + W"ti) ■ ^ ' 

Note also that this result for Pi is not the accurate stationary solution, since the fact 
that it does not vanish for large r invalidate the truncation of the space to the the first 
two phonon states. To find the true stationary solution we have to include more states, 
this can be easily done numerically. The result will be discussed in the next Section, 
but it is clear that the system cannot be cooled at very low temperature in region D. 
It is also clear that going even more deeply inside the sequential tunneling regime the 
situation will not improve, and finally one would find that electronic transport heats 
up the oscillator. 



4 Numerical results 

In order to confirm the analytical results and to investigate their limitations in this Sec- 
tion we discuss the numerical solution of the full master equation ((5]). This can be done 
by truncating the dimension of the space to 2pjvf states such that the eV ,eVg,kBTm, 
kgTe <^ pmTT'^- In practice in the following pm = 30 will be sufficient for our purposes. 

We begin by considering an ideal case of very cold electrons ksTe <^ Tilu, and 
long relaxation rate of the phonons tFq/H — 100, with Fo/hij = 0.1 and A = 1. The 
resulting effective temperature of the oscillator is shown as a function of the bias and 
gate voltage in Fig. [S] The plot shows how important are the bias condition, since 
the effective temperature dependence clearly allows to recognize the different regions 
we discussed in Section [S] since the temperature changes abruptly at their borders. 
The results confirm the analysis of the previous Section, regions A, B, and the stripe 
of C regions correspond to a stationary solution with very low effective temperature 
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Fig. 3 (Color online) Effective temperature of the oscillator in the eVg-eV plane. In this 
calculation Te/huj = 0.01, Tm = fiu), t/Iq = 100, A = 1, and Fo/liJ = 0.1. We assume 
symmetric junctions. In this figure it is easy to recognize the different regions shown in Fig. |2] 
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Fig. 4 Effective temperature as a function of the coupling constant A for the same parameters 
of Fig. [3] for different bias points. Specifically the full, dotted, and dashed lines correspond to 
the bias points {eVg,eV)/hu> = (0,0.5), (0,2.5), and (1,0.5), respectively. 



(cfr. Fig. [3]). The effective temperature is much smaller than the mechanical reservoir 
temperature. In particular the lowest temperature is obtained for any of the C regions, 
regardless of the value of the bias voltage. We also can see that the temperature is 
higher in the D rhombi for the reasons explained in Section [3] 

In Fig. |4]we plot the coupling dependence of the effective temperature for the same 
parameters of the Fig. [3] at different bias points as indicated in the caption. The full 
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Fig. 5 (Color online) Effective temperature for Te/hui = 0.1, Tmlfiijj = 5, t/Fq = 100, A = 1, 
and Fo/u! = 0.1. Compared to Fig.[3]we plot here a much larger region of the plane, since the 
temperature is higher. 



and dotted line correspond to the case A, and C, respectively. We see that the optimal 
cooling is attained for A = 1, as found in the analytic calculation. The dashed line 
shows the behavior in the middle of the Coulomb diamond, and for bias voltage larger 
that 2hij, that constitutes the inelastic threshold in this case. We see that an increase 
of the coupling leads to a strong heating of the device, which is expected since the 
electrons have enough energy to release one phonon in the oscillator every time they 
traverse the device. For very large coupling the effect is suppressed by the reduction of 
all the rates due to the exponential vanishing of the overlap of the harmonic oscillator 
wavefunctions. 

Finally let us consider the case of Tm = 5hijj and still cold electrons Te = O.lhui. 
The effective temperature for this case (the other parameters are the same of Fig. |3} is 
shown in Fig. (5] The main feature is clearly the presence of a deep trench at the border 
of the conducting diamond. Again this is in agrement with the analytical results, that 
shows that Te// should not depend on eV along that line. We also see that in the 
conducting region for large enough voltages the current through the device heats up 
the oscillator, as is usually the case. On the other side, too far from the conducting 
region the temperature is not modified since transport is completely suppressed. 

As a final check we also verified that when the electronic temperature Te coin- 
cides with Tm, cooling is not possible. There is always a dependence of the effective 
temperature of the stationary state on the bias condition, and in some case we find a 
very small reduction of the temperature of the oscillator, but clearly in this regime the 
device could not be used as a cooler for the oscillator. 
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5 Conclusions 

We have shown that a molecular (single electronic level) single-electron transistor cou- 
pled to an oscillator can act as a cooler for the oscillator if the electronic temperature 
is lower than the temperature of the environment to which the oscillator is coupled. 
The heat transport mechanism is not trivial since it involves the phonon absorption 
and emission, that can be finely tuned by the bias condition of the single electron tran- 
sistor. We find that optimal cooling is achieved by biasing the device near the border 
of the sequential tunneling conducting region of the Vg-V plane. By using for instance 
normal-metal-superconductor tunnel junctions electronic cooler, it is possible to gen- 
erate cold sources of electrons and then to exploit them to cool nano-electromechanical 
oscillators. 
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